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ABSTRACT 

Our new formulation of the Spherical Collapse Model (SCM-L) takes into account the pres- 
ence of angular momentum associated with the motion of galaxy groups infalling towards 
the centre of galaxy clusters. The angular momentum is responsible for an additional term 
in the dynamical equation which is useful to describe the evolution of the clusters in the 
non-equilibrium region which is investigated in the present paper. Our SCM-L can be used 
to predict the profiles of several strategic dynamical quantities as the radial and tangential 
velocities of member galaxies, and the total cluster mass. 

A good understanding of the non-equilibrium region is important since it is the natural 
scenario where to study the infall in galaxy clusters and the accretion phenomena present in 
these objects. Our results corroborate previous estimates and are in very good agreement with 
the analysis of recent observations and of simulated clusters. 
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1 INTRODUCTION 

The cosmic structures we observe in the Universe are thought to 
be originated by small perturbations in the primordial matter dis- 
tribution, enhanced in time by self-gravity. Density perturbations 
decouple from the Hubble flow and drive the collapse of matter 
into structures (groups and clusters of galaxies) through a highly 
non-linear process. In a Universe with cold dark matter and a cos- 
mological constant (A-CDM), structures of smaller scale decouple 
before structures of larger scales do: this means that galaxy clusters 
are built up through merging and coalescence of smaller galaxy 
groups. This "bottom-up" scenario has been widely tested using 
large A'^-body hydrodinamical simulations and is consistent with 
the most recent measurement of the cosmological parameters. 

An accurate general description of the cosmic structure col- 
lapse is analitically inpraticable. To tackle the problem, simpli- 
fying assumptions are required. The most common and accepted 
assumption is to consider only spherically symmetric density per- 
turbations: in this case, each perturbation can be described as an 
independent Friedmann universes on its own, with constant posi- 
tive curvature, which first expands and then collapses mantaining 
its spherical shape. This Spherical Collapse Model (SCM) was first 
developed by |Gunn & Gott| ( (T972) , [Slik| ( (T974l ), and [Gurm| ( [T977l >, 
and it has since became a widely accepted approach in the liter- 
ature (e.g. |Peebles|[T976l [T980l [Schechte r 1980; S hectman||1982| 
[Ostriker et al.|1998^ . Many authors demonstrated that despite its 
simplistic assumption the SCM provides a quite good description 
of the dynamics in the outskirts of clusters, where the matter is not 
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in energy equilibrium and is characterized by an overall infall mo- 
tion towards the cluster centre (e.g. 'Prada et al.'2006' 'Betancort- 
|Rijo et al.f "2Q06 , Bilic ki & Chod orowski 2008 , Cupani, Mezzetd] 
& Mardirossian 2008 , hereafter CMM08). Corrections to the SCM 
are generally required in the cluster cores, due to the interpenetra- 
tion of different collapsing shells of matter (an effect referred to as 
shell crossing; see e.g. |Sanchez-Conde et al.|2007| and references 
therein). 

One of the major shortcomings of the SCM is that in princi- 
ple a purely spherically symmetric scenario does not allow for any 
non-radial motion of the collapsing matter. In fact, galaxy dynamics 
in observed clusters show a significant amount of tangential motion 
( |Ryden & Gunn|1987p . [Diaferio & Geller|(T997) demonstrated that 
the amplitude of the caustic surfaces observed in the redshift-space 
distribution of galaxies is related to the galaxy escape velocities, 
which are by definition non-radial (see e.g.lDiaferio 1999 , R ines &| 
|al.|2000||20UT| |2003 ; Rines & Diaferio 2006; and Diaferio 2009 for 
a review). More recently, _Cupani, Mezzetti, & Mardirossian_^20l"0^ 
[hereafter CMMIO} showed that radial velocities and tangential ve- 
locities of galaxies are statistically of the same order of magnitude 
in modulus. Observation of line-of-sight velocities can therefore be 
used to estimate the infall pattern and the mass profile of clusters 
up to the outermost regions (corresponding to radii exceeding seven 
times the virialization radius). 

In this paper we take into account the effect of tangential mo- 
tion into the SCM. Considerations based on the A-CDM bottom- 
up scenario suggest that angular momentum in cluster outskirts is 
mainly associated with galaxy groups which belongs to the clus- 
ters, rather than single galaxies. In fact, galaxies rarely infall iso- 
lated, since groups appear to be formed long before they merge 
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into clusters. We will show that the tangential motion of infalling 
galaxies is not only due to the intrinsic velocity dispersion of the 
infalling groups which the galaxies belong to, but also to the cou- 
pled effects of the angular momentum conservation and the torque 
produced by the presence of nearby structures. Several authors dis- 
cussed the effect of the angular momentum in cosmic structures, 
but they mainly focused on the growth of the total angular momen- 
tum and/or on the behaviour of the inner, relaxed regions (e.g.lRy 
den & Gunn||1987[|Ryden|[l988][Cateian & Theuns][l996l [Aviia 
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Reese et al.|1998[|Nusser|2001[|Ascasibar et al.|2004[|Del Popoio 
& Kroupa||2009 i. On the contrary, we will not take into account 
here the overall angular momentum of the galaxy cluster which can 
arise from the misalignment between the inertia tensor and the tidal 
forces acting on the collapsing structures i White|1984p , nor the pro- 
cess of violent relaxation ( [Lynden Bell| 1967[ >, which transforms 
the infall energy into non-radial velocity dispersion on a collapse 
timescale. But we will follow a different approach and focus on the 
non-equilibrium region of clusters, where the presence of angular 
momentum is expected to slow down the overall infall motion of 
matter We will show that the SCM with angular momentum (here- 
after SCM-L) better describes, compared to the standard SCM, the 
infall velocity profile and the mass profile in the outskirts of clus- 
ter; in fact it predicts shallower infall velocity profiles which are 
consistent with the results of CMM08 about the velocity profiles. 

We wish to stress that considering the effects of the angular 
momentum in the dynamics of the outskirts of clusters necessarily 
implies an approach which is linked to the dynamics of the separate 
structures infalling towards the cluster core, so we have not to deal 
with the cluster total angular momentum. 

In the following, we will discuss how to to take into account 
the effect of the angular momentum and how to derive it from 
the power spectrum of the density perturbations which originated 
the clusters (sectionlSj. The derived equation provides information 
about both the matter distribution and the kinematical properties 
of clusters. We will thoroughly compare these predictions with our 
results pubhshed in CMM08 and CMMIO (section|4]l and discuss 
their implications in cosmology (sectionBl. The details of the cal- 
culation are generally omitted in the main text, and a full account 
on them can be found in the appendices. 

In the present paper we adopt Ha = 70 km s^^ Mpc^^. 



2 INSERTING ANGULAR MOMENTUM 

In the following, we detail the SCM-L scenario. We consider a 
spherical region with centre O, surrounded by an (almost) isotropic 
Universe with its own fluctuations. We assume at first that this re- 
gion is actually completely homogeneous, i.e. there is not a cluster 
in it. Let us consider an object (within a mass shell) at a distance a 
from O, with a tangential velocity v_i_. The evolution of this v±_ is 
described by the evolution of the specific (i.e. per unit mass) angu- 
lar momentum C with respect to O; 
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where aj-'x is the specific torque and J^^ is the specific transversal 
force. 

Let us suppose now that in the same spherical region there 
is a galaxy cluster centered in O. In this case the mass shell, that 
was at a distance a from the centre in the previous case, is now 
at a distance r from O. The specific angular momentum C of an 
object in the mass shell at the distance r from the center and with 
transversal velocity v± will evolve according to 



l{avi_) = aJ^x, 



(1) 



where v±, C, and .Fx are barred to indicate that in principle they 
may be different with respect to the previous case. 

The velocity of an infalling system is perturbed by the force 
produced by external structures. Structures on scales comparable 
to, or larger than the cluster scale will affect the overall motion of 
the cluster and of its outskirts, in such a way that the relative ve- 
locity between the cluster and its outskirts will not change. On the 
other hand, structures on scales comparable to, or smaller than, the 
group scale are not big enough to significantly affect the infall mo- 
tion of groups. So, only structures in a range of scales will affect 
the relative velocities of the infalling systems. As we will later jus- 
tify, the dominant contribution comes from scales on the order of 
~ 10/i~^ Mpc, and the difference between r and a are of the order 
of a Mpc, or less. We will therefore neglect the difference between 
J^x and T±_ inside the cluster, and we will assume they are equal, 
making the problems analytically tractable: 



(3) 
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In order to estimate C = av±^ in the no-cluster case, we can 
use the peculiar velocity v, obtained through the continuity equa- 
tion linking the density contrast 5 \.ov (see, e.g., |PeacockH199^ ). 
The ortogonal velocity v±_ is obtained as the projection of v onto a 
plane perpendicular to the radial direction. Under the isotropy as- 
sumption, this gives v±_ = ^2{v^)/3, while radially this gives a 
contribution |?;[| | = ^{v^)/3 to the infall velocity. 

The knowledge of £ = rv± allows to insert it into the Fried- 
man equation as follows: 

(4) 
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where p is the mean density of the perturbation within r, G is the 
gravitational constant, A is the cosmological constant, and c is the 
speed of light. Equation (HI gives an exhaustive description of the 
SCM-L. 



3 PROCEDURE 

Three main statuses can be recognized in the evolution of a density 
perturbation: 

(i) A primordial status, when the pertubation uncouples from 
the Hubble flow and expands at a slower rate than that of the un- 
perturbed Universe; 

(ii) A turnaround status, when the perturbation reaches its max- 
imum amplitude and starts collapsing; 

(iii) A virialization status, when the collapse eventually ends 
and the dynamical equilibrium is reached. 

These statuses are reached at different times, and they depend on 
the primordial overdensity of the perturbation which originated 
the cluster. Overdensities which reach either the turnaround sta- 
tus or the virialization status at the end of their evolution are par- 
ticularly interesting for studying the evolution of the clusters (see 
e.g. [Navarro, Frenk, & White T995; 1996" 1997, Bulloc k et aT 
|2001||Lokas & Mamon 2001 ; Cupani, Mezzetti, & Mardirossian 



2008, 201 0, hereafter CMMIO I. The SCM allows us to analyze 
these two dynamical statuses stopping the cluster evolution at dif- 
ferent times (which can or cannot coincide with the present time 
to). We use throughout the subscripts t and v to refer to the 
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turnaround stage and to the virialization stage, and the subscripts 
i, f , and to refer to the initial time and the final time of the evolu- 
tion, and to the present time, respectively. 

The cosmology is defined in agreement with the A-CDM pre- 
scription: 



Q.M = ^h{z) + Q,r,M{z)\ VlM{z)+Q.f^{z) = 1, 



yz: 



flM{z) 



»M.0(1 + Z)^ 
f2A,o + ^M,o(l + z)^ 



(5) 

(6) 



Here f^x denote the density fraction of the species X with respect 
to the critical density. The subscripts are defined as follows: M = 
matter, b = baryonic matter, DM = dark matter, A = dark energy 
or cosmological constant. We adopt throughout the concordance 
value for the baryonic density fraction, flb,o = 0.04 (e.g., |Borgani| 
|et al.|1997l|Borgani et al.|2004| , and we also choose many other 
values for the matter density fraction, namely Qm,o ~ 0.20, 0.25, 
0.30, 0.35, and 0.40. We use SIm.o values usually rejected by the 
literature, just in order to investigate the robustness of our results 
for different nM,o- 

The power spectrum of the perturbations P{k) (with k the per- 
turbation wave number) is defined, according to the prescription of 
[Eisenstein & Hu| ( |1998| >, using their transfer function T{k) (accu- 
rate to better than 5 per cent for f2b,o/f2M,o ^ 0.5): 
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n{k)- P{k) = kT{kY. (7) 



Here Tb(t) is the component associated to the baryonic matter and 
7bM(fc) is the component associated to the dark matter. The re- 
sulting power spectrum is normalized to the root-mean-square den- 
sity fluctuations within a sphere with radius rg = 8/i~^ Mpc, as, 
imposing, as in the reference simulation of [Borgani et al.|[2004| 
as, = 0.8, in agreement with the value obtained by the latest mea- 
surements of the Wilkinson Microwave Anisotropy Probe (WMAP) 
(see e.g. |Komatsu et al.|2008| ): 
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Both the initial overdensity 5i and the angular momentum £ 
of the perturbation are extracted from the power spectrum P{k). 
To compute 5i (within ri), we apply a gaussian filter to the power 
spectrum, with radius tq = l/i~^Mpc, since we are interested on 
progenitors at the galaxy cluster scale (see appendix [All. To com- 
pute £, we use a band-pass filter, peaked at a scale of ~ lO/i^^ 
Mpc, to isolate the contribution at the scale of galaxy groups (see 
appendixlBl. The last band pass filter is designed on the assumption 
that galaxy groups are mainly responsible for the angular momen- 
tum term which has to be introduced in the standard SCM. This as- 
sumption is consistent with the widely accepted bottom-up scenario 
(which looks at galaxy groups as the building blocks of galaxy clus- 
ters), and it will be justified a posteriori by the dynamical analysis 
of the catalogue of simulated galaxy clusters we use (see section 
|4.2[ (. The use of both the gaussian filter, both the band-pass filter, 
cuts away almost all the non-linear part of the power spectrum. For 
this reason, in the present paper, we used the linear P{k). 

The evolution of the perturbation is described by time integra- 
tion of equations dSb and iH. We need as input 5i and £,{ in order to 
obtain as output R{ = rt/ai, d_R/dr|j, Af = (1 + SiJR^'', and 
£.{. The integration is performed between an initial time ri and the 
final time T{, corresponding to initial redshift value Zi = 1000 and 
the final value zt, which depends on our choice of the end of the 



evolution, i.e. Z[ = 0, 0.1, 0.2, 0.3, and 0.5, which correspond to 
possible redshifts of observed clusters, r is defined as 
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according to the adopted cosmology. 

We focused on the cases of a perturbation reaching either the 
turnaround state or the virialization state at the end of its evolution. 
The turnaround state is defined by the condition di?/dr|j. — 0; 
we used this condition to obtain the initial and final radius, Vi^ and 
rf,t, and the initial and final overdensity, 5i^t and (5f,t, of a perturba- 
tion evolving into the turnaround status just at T{. The virialization 
status is quite ambiguous, since neither the SCM nor the SCM- 
L are able to describe the details of final virialization, when shell 
crossing occurs. However, it is usually assumed in the literature 
that virialization occurs when, in the SCM scenario, the overdensity 
collapses into a point. In the SCM-L scenario, angular momentum 
prevents the collapse, and we have assumed for analogy that virial- 
ization corresponds to the state when Rf reaches its first minimum. 
We call ri,v and (5i,v the initial radius and the initial overdensity of 
a perturbation which corresponds to this condition. 

We computed the turnaround radius, Rt = rf.t/af, and the 
turnaround overdensity. At = l + 5f,t, of this evolved perturbation. 
The ratio between turnaround mass and virialization mass is 



m ^ 1 + 5i,t 

Afv ^ 1 + (5i,v 



(10) 



while the ratio between turnaround radius and virialization radius 

is 



n ^ f Ay A/t 
rv ~ V ^t,t Mv 



1/3 



(11) 



where Av = [IStt^ + 82(»m,o - 1) - 39(nM,o - lf]/^Mfi is the 
standard virialization overdensity (e.g., [Bryan & Norman|1998^ . 

We follow the evolution of an overdensity profile from ri to rf . 
We mapped the initial overdensity profile as a set of pairs {ri, 5i}fc, 
with A: = 1, . . . , 10, each one describing a single perturbation, and 
let the pairs evolve into {Rt, Af}fc, to obtain the final overdensity 
profile. We are also able to compute the final profile of the radial ve- 
locity ii|| f and the tangential velocity ii_L,f (see appendixlBl. In the 
analysis, we focused on the following quantities, which are relevant 
for the dynamics on the non-equilibrium region of galaxy clusters 
(CMMIO): 
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where the tilde denotes normalization to turnaround values, as sug- 
gested by CMM08. 

The results of our analysis are discussed in the next section. To 
assess the difference between the SCM-L and the standard SCM, 
we switched on and off the angular momentum term in the evolu- 
tion equations ll3l and (HI. 
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Table 1. Ratios between the SCM and the SCM-L outcomes for Rt 
^t,SCM/^t,SCM-L. for different choices of SIm.o and Zf. 
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Table 2. Ratios between the SCM and the SCM-L outcomes for At, i.e. 
^t.SClvl/At, SCM-L for different choices of r2M,o and Zf . 
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RESULTS 













In section 14.11 we show that the SCM-L better estimates the 
turnaround overdensity, while no large differences are seen in the 
SCM and SCM-L ratios between the turnaround radius and the viri- 
alization radius, and in the SCM and SCM-L ratios between the 
turnaround mass and the virialization mass. 

In section [42I we show that both SCM and SCM-L give quite 
good results in the determination of the cluster mass profiles, but 
only the SCM-L succeeds in a good fit of the radial velocity profile 
of the cluster components. 



4.1 Turnaround 

We first discuss the predictions of our model at turnaround. The 
ratios between the different values of Rt and At obtained from 
the SCM and from the SCM-L are listed in table [T] and in table 
[2] The SCM-L predicts smaller turnaround radii consistently with 
higher turnaround overdensities. This means that the presence of 
angular momentum slows down the infall motion of cluster com- 
ponents in its outskirts. In fact, a lower infall velocity implies that 
the turnaround region has to be shifted inwards, to higher-density 
regions. This effect is remarkable and is consistent with the results 
of CMM08 about the velocity profiles, obtained from the statisti- 
cal analysis of simulated clusters. The ratios between the SCM and 
the SCM-L outcomes for Rt and At, i.e. -Rt,scM/-Rt,scM-L and 
At, scm/ At, SCM-L, are quite similar for different choices of f^M.o 
and Zf (see tables[T]and[2]l. 

In order to show the dependences of the obtained values of 
At, Mt/Mv, and rt/r\ on SIm.o and Zf, we obtained the following 
analytical fittings, which reproduce the predictions of the models 
with an accuracy better than 3 per cent for 0.2 ^ S7m,o ^ 0.4 and 
sC Zf sC 0.5: 
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Values obtained from equations ([151, l|T6j, and iTTtJ in the above- 
mentioned ranges of Q,Mfi and Zf are represented with a grey scale 
in figure (??) for both the SCM and the SCM-L (upper and lower 
panel, respectively). The difference between the two models mostly 
affects the values of At, while it is less evident for the ratios 
J\ft/A/v|fit and rt/rv|fit. This indicates that the turnaround scale 
and the virialization scale are similarly affected by the introduc- 
tion of the angular momentum term into the evolution equation. It 
means that the values of Mt/M^ and rt/r^ obtained by CMM08 
in the SCM scenario do not differ significantly from the values one 
should obtain in the SCM-L scenario. 

Both SCM and SCM-L are in quite good agreement with 
the observations and the analysis of cosmological simulations for 
f^M.o = 0.3 and zt = 0. CMM08 evaluated a turnaround over- 
density At = I6I3, taking into account cosmological simula- 
tions porgani et al.|1997| r2004 , Biv iano et al.|2006^ , which was 
only marginally consistent with the SCM scenario. This discrep- 
ancy fosters the necessity of considering the angular momentum in 
the SCM. Our new SCM-L succeeds to be in quite good agreement 
also with |Rines & DiaferiO| (2006l (who found a ratio between the 
turnaround mass and the virialization mass equal to about twcPjl. 



4.2 Profiles in the non-equilibrium region 

We now consider the mass and velocity profiles predicted by of the 
SCM-L with JIm.o = 0.3 and zi = 0, coherent with the concor- 
dance A-CDM cosmology and the near cluster catalogues. To bet- 
ter appreciate our results, we compare our theoretical profiles with 
those obtained from a catalogue of 9631 member galaxies iden- 
tified within 1 14 simulated clusters by [Biviano et al.| P006| ) (see 
[Borgani et al.|199T] and [Borgani et al.|2004 for details on the cos- 
mological simuation; the catalogue is described also in CMM08 
andCMMlO). 

The profiles of M — M/Mt versus R = R/Rt are shown 
in figurefTl The theoretical profiles predicted by SCM (bold dashed 
line) and by the SCM-L (bold solid line) are superimposed to the 
mean profile obtained from the distribution of the 9631 member 
galaxies (grey solid line); the la region of the distribution of cluster 
mass profiles is also shown (shaded region). Both the SCM and the 
SCM-L are in very good agreement with the mean simulated pro- 
file. Nevertheless, the SCM-L turns out to be slightly more accurate 
to describe the mass distribution beyond the turnaround region. 

The NFW mass profile by .Navarro, Frenk, & White] jl995| 
|1996[ |1997^ is also plotted in figure [2] (narrow dashed line). As 
already pointed out by CMM08, the NFW profile seems not to suc- 
cessfully describe the non-equilibrium region of galaxy cluster. We 
take into account also the modification of the standard NFW profile 
suggested by|Tavio & al.|p008} (narrow solid line). Evidently, this 



^ The original result is A/t/Af200 = 2.18ib0.19, where M200 is themass 
enclosed in a sphere with overdensity A = 200. Assuming that Mv ~ 
l.lAf200, according to what stated by |Rines & Diaferio, 1 ,2006) , we obtain 
Mt/Mv = 2.0 ±0.2. 
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Figure 1. Comparison between mass profiles. Bold dashed line: SCM pre- 
diction; bold solid line: SCM-L prediction; narrow dashed line: NFW pro- 
file ^Navarro, Frenk, & W hitell995| |1996][T997,) ; narrow solid line: modi- 
fied NFW profile (Tavio & al. 2008). The grey solid line and the shaded area 
represent the mean profile and the \a region as extracted from a catalogue 
of 9631 member galaxies ^Borgani et al.|1997l|Biviano et al.|200'6). 



Figure 2. Comparison between radial velocity profiles. Bold dashed line: 
SCM prediction; bold solid line: SCM-L prediction; narrow dashed line: 
Yahil approximation jYahil| 1985,1 ; narrow solid line: Meiksin approxima- 
tion jVillumsen & Davis|1986} TThe grey solid line and the shaded area 
represent the mean profile and the lo" region as extracted from a catalogue 
of 9631 member galaxies ^Borgani et al.|1997||2004[|Biviano et al.|2006) . 



profile is in better agreement with the simulated cluster distribution 
and is fairly consistent with our SCM-L prediction; it evidences 
a significant convergence between [Tavio & aTlpOOS^ and our ap- 
proach. 

The radial velocity profiles of the cluster components are 
shown in figure [2] Once again, the theoretical profiles predicted 
by SCM (bold dashed line) and by the SCM-L (bold solid fine) are 
superimposed to the mean profile obtained from the distribution of 
the 9631 member galaxies (grey solid line); the Icr region of the 
distribution of cluster mass profiles is also shown (shaded region). 
Here the difference between the SCM and the SCM-L stands out 
clearly, in fact the two models predict different values of Af at 
turnaround. In particular, the SCM gives lower values of At than 
the SCM-L does, with a consequent steeper velocity profile in the 
SCM scenario. The SCM systematically overestimates the amount 
of radial velocity in the whole overdensity range from the extreme 
outskirts of clusters to the virialization core. This is not the case of 
the SCM-L scenario, due to the presence of the angular momentum 
term, which has been shown to be responsible for slowing down the 
infall motion in the non-equilibrium regions. The abovementioned 
agreement between the SCM-L prediction and the simulated mean 
profile gives a strong support to the hypothesis that angular momen- 
tum, due to the motion of galaxy groups in cluster outskirts, has to 
be taken into account. The agreement also justifies a posteriori our 
assumption that the dominant contribution to angular momentum 
comes from scales of the order of ~10/i~^ Mpc (see sectionl2|. 

Two different literature approximations for the radial velocity 
profile are also shown in figure [2] namely the Yahil approximation 
Vy (|Yahil|1985[ narrow dashed line), 



and the Meiksin approximation I'm jVillumsen & Davis|1986[ nar- 
row solid line). 



i)M(A) 



\/3 A-1 

3 (A + 2)1/2 



(19) 



wy(A) = 



3A1/4 



(18) 



The Yahil approximation is in good agreement with the profile pre- 
dicted by the standard SCM for A < 30, while the Meiksin approx- 
imation is nearly coincident with our new profile predicted by the 
SCM-L, and thus it seems to be the best choice. It confirms what 
was found by CMM08. 

We also computed the ratio between the tangential velocity 
and the radial velocity predicted by the SCM-L, Kv = u_l/ii||, as 
a function of the radial coordinate R = R/Rt (where Rt is the 
turnaround radius). The profile obtained from this latter function is 
shown in figure [3] superimposed to the mean profile obtained from 
the distribution of the 9631 member galaxies (grey solid line); the 
1(7 region of the distribution is also shown (shaded region). The 
quantity K^j is evidenced by CMMIO to be crucial to relate the 
mean observed line-of-sight velocity of member galaxies in clus- 
ters to their mean infall velocity (towards the cluster centre). The 
knowledge of the mean infall velocity profile is crucial to estimate 
the mass profile (see CMMIO for the details). 

The SCM-L prediction of Kv agrees within the Icr uncertainty 
with the mean profile traced by the simulated galaxies, and it is 
nearly constant along the radial direction, where the value of Kv 
is about 0.95 in the turnaround region. Kv = 0.95 is the fiducial 
value adopted by CMMIO to identify the 'fair region' in the redshift 
distribution of galaxies, i.e. the region where the galaxies are shown 
to better trace the infall motion, and subsequently to better obtain 
the cluster mass profile. 



6 Guido Cupani, Marino Mezzetti, Fabio Mardirossian 



^ 1: 




Figure 3. Ratio between tangential and radial velocity. Bold solid line: 
SCM-L prediction. The grey solid line and the shaded area represent the 
mean profile and the Icr region as extracted from a catalogue of 9631 mem- 
ber galaxies ^Borgani et al.|1997l|Biviano et al.|2006). 



5 CONCLUSIONS 

Our new formulation of the Spherical Collapse Model (SCM-L) 
takes into account the presence of angular momentum associated 
with the motion of galaxy groups infalling towards the centre of 
galaxy clusters. The angular momentum (estimated from the power 
spectrum of density perturbations) is responsible for an additional 
term in the dynamical equation which is useful to describe the evo- 
lution of the clusters in the non-equilibrium region. Once the evolu- 
tion equation is solved, our SCM-L provides a complete description 
of the dynamical state of the perturbations at any given redshift z 
as a function of the adopted values of matter density function Slo- 

We run our model with different values fio between 0.2 and 
0.4 and redshift z between and 0.5. In each case, we computed 
the turnaround overdensity At and the ratio between the turnaround 
scale and the virialization scale. 

Our SCM-L can be also used to predict the profiles of sev- 
eral strategic dynamical quantities (radial and tangential velocities 
of member galaxies, and total cluster mass). The agreement be- 
tween our SCM-L and the simulations is always remarkable in the 
whole non-equilibrium region (spanning approximately from 0.5 to 
2 turnaround radii rt), which is the region of interest of the present 
paper. 

The non-equilibrium region is the natural scenario where to 
study the infall in galaxy clusters and the accretion phenomena 
present in these objects. The method described in the present pa- 
per can be also generalized to alternative cosmological models. 

Our results corroborate our previous estimates (CMM08); 
moreover, they are in very good agreement with the analyses of 
recent observations ( [Diaferio & Geller|1997^ and of simulated clus- 
ters (see, e.g., CMMIO). 



, Hernandez X., 1998, ApJ, 505, 37 
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APPENDIX A: PRIMORDIAL OVERDENSITY OF THE 
PERTURBATIONS 

This appendix provides the recipe for computing the initial over- 
density profile of a spherically symmetric density perturbation from 
the power spectrum. We use the results of Bardeen et al.| ( |1986[ ), 
Lahav et ar](|1991), and|Lilje & Lahav| ( [T991| ); see also |Cupani,| 
Mezzetti, & Mardirossian J2008 1 for an example of the same ap- 
proach. 

The overdensity profile Si (within ri) is obtained from the 
statistic of peaks in the density field described by the power spec- 
trum. To identify the progenitors at the cluster scale, we apply a 
gaussian filter with radius re ~ l/i~^Mpc to isolate the contribu- 
tion of the galaxy cluster scale: 



PG(fc) = P(fc)e'(''G'=)V2. ^^jj 

Pg (fc) is represented as a dashed line in figurep] the original power 



spectrum P{k) is also represented for comparison (dotted line). We 
denote with af the {I + l)-th even moment of Pg(A;): 



cri = 



27r2 



PG(fc)fc2('+iMfc, 



(A2) 



If i/ao is the expected height of a peak, we can express the differ- 
ential number density of peaks as a function of v. 



N{v 



(27r)2r 






D ' —1 



H(7,7!/), 



ol 



CT0O"2 Cr2 



H(7,w) 
A(7) 



TO^ + 3-y^w + A{'y)w^ + B{-y) e 



-C{-,)n 



1 + D(7)e-^(T)'' 



432 



/l07(9 - 572)5/2 ' 

5(7) = 1.84 + 1.13(1 - 7^)^■''^ 

5 
^'''^^ ^ 2(9-572)' 

D(7) = 8.91 + 1.27e^'5'^', 

£(7) = 2.58e'-''^^'. 



(A3) 
(A4) 

, (A5) 
(A6) 

(A7) 
(A8) 

(A9) 
(AlO) 



The linear approximation for the present-day overdensity profile is 
consequently obtained as 

3 



5lm,o(r) = „ 2 / 



kji{kr)Pf,i{k) X 



u — -y^u — jO Or'i 2 

^37(1-7^) 



1-72 



dk. 



/- e(^)AA(^)d^ 



(All) 
(A 12) 



e(^) = 3(l"7') + (l-216-0.97^)e-<^/^»^-/^>\ ^^^3^ 
\/3(l-7^) + 0.45+(^)V^ 

We assume that at the initial time the perturbation is not yet 
uncoupled from the Hubble flow, and that the initial overdensity 
is thus well predicted by the linear approximation: 5i = Sn^,!- 
This requirement is generally satisfied at high redshift: Zi ~ 1000. 
To compute the initial profile 5iin,i(r), we simply revert back the 
evolved profile (5k„,o (*") using the growing-mode linear solution of 
the SCM: 



5lin,i{r) = 5\in,o(r) 



DAZr 



D+{0) ' 



(A14) 



D+iz) 



5Qm{z) 
2{l + z) 



Qm{z)^^'' - nA(z) + 



+ri+^)(i+"^ 



(A15) 



APPENDIX B: FRIEDMAN EQUATION WITH ANGULAR 
MOMENTUM 

This appendix describe the technical details concerning the evolu- 
tion of a spherically symmetric density perturbation with angular 
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momentum. We first discuss how to obtain the angular momentum 
from the power spectrum of perturbations and then detail the rou- 
tine to numerically solve the Friedman equation |4l. 

The angular momentum can be computed from the peculiar 
velocity field obtained from the power spectrum. To isolate the an- 
gular momentum contribution of galaxy groups, we adopt the fol- 
lowing band-pass filter, using a scale radius r^ for the cluster scale 
and a scale radius r„: 



Pi,p{k) = P{k) 
f 3 



1- 



sin(fcrc 



kr-ci 



sin(fcr 



fcrgi- cos(fcrg 



(kr^rr 



(Bl) 



Phpik) is represented as a dashed line in figure Hi the original 
power spectrum P{k) is also represented for comparison (dotted 
line). The radii r^ and r^^ are computed as the mean separation 
between neighbouring clusters and neighbouring groups, respec- 
tively, defined as functions of the cluster number density Uci and of 
the group number density n^,.: 



r-ci 



4-Knr 



1/3 



4-Kn„ 



1/3 



(B2) 



The cluster number density is extracted from the large cosmologi- 
cal simulation of Borgani et aLl j UggVlpOO? *) (see also |Biviano et al.| 
|2006; CMM08; CMMIO) . The estimate for galaxy groups is taken 
from |Girardi & Giuricin| l |2000) , and corresponds to the number 
density of objects with mass larger than 9 x lO^^h^^ Mq (...) 



nci~ 1.6 X 10^^/i^Mpc"^ 



rici ~ 1.0 X 10"^/i^Mpc' 



rci = 24.6/i"^ Mpc; (B3) 

^ rci = 6.20/i"^ Mpc. (B4) 



The peculiar velocity field is obtained through the continuity rela- 
tion between S and v. If the motion of the group is only due to large 
scale perturbations, S is well predicted by the linear approximation. 
We obtain in this case 



{v') = 



fs 



27r2 

(1+Z)2 

FiS) _ 



Php{k)dk 

" k^ P^,p{k) al dk _ 
27r2 FT' 





dlnD+ 
dlnao 



(B5) 



(B6) 



The ortogonal velocity v± is obtained as the projection of v onto 
a plane (cf. section [2](. Under the isotropy assumption, this gives 
v± — \j2{v'^) /Z. From equation i B5 i we finally obtain 



C(z) 



av± — HoaoKQ{z) 



K = 



Q{z) = 



I2 roo 


fc3Pbp(fc) al 


dfe 


27r2 fc2 


k 
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(B7) 



(B8) 



(B9) 



The knowledge of C allows the estimation of C from equation 
1I3J. Equation (HI and equation ll3j have to be solved together, since 
both the radius r and the angular momentum C change in time. 
After rearranging, we obtain 



R{r) = -244^M + [(1 - «M.o) - ^ 



R{t) + 



«M,o[l + (5i(n)] 
2i?(r)2 

£{t) = KRQ{t); 
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T = Hot^ 
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; (BIO) 
(BID 

-3/2 

(B12) 
(B13) 



where the dots denote first- and second-order derivatives with re- 
spect to r. The boundary counditions needed to solve the system 
| |B10| >-i |bT7) are set by assuming that at the initial time the pertur- 
bation is not yet uncoupled from the Hubble flow, i.e. r^ = a^, 
but the peculiar radial velocity is given by the linear reltion. This 
requirement is satisfied at high redshift: zi ~ lOOO. Consequently, 



Ri 



1/2, 



Ci = KQi. 



n3/2 



(B14) 
(B15) 
(B16) 



Once solved, the system i |B10[ i-( [BTT) provides the values of 
R{, R{ and C{ for any given value of Tf > n. All the relevant 
physical properties of the perturbation, including the infall velocity 
Vr,f, are consequently obtained as follows (r^ is given in comoving 
units): 



ri = Rin; 

5f = (1 + 5i)i?r^ - 

"II. f ^ yrHoRfri-, 
Hi 
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